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Int roduce .on and S u n-nary 


The final report documents the research work accomplished on the NASA 
Crant NSC 8037 entitled "Atmospheric Effects oa Coherent Laser Systems" 
during the period June 1978 through May 1979. This report includes the 
material from the semi-anr.ual report along with the acditional analysis 
developed during the latter half of the Crant period. Most of the material 
from eq. (16) through eq. (33) is due to the research since December 1978. 

The objective of the research is to investigate the effects of che 
truncation of the gaussian beam by the telescope on the signal-to-noise 
ratio of the NASA Laser Doppler System. The analysis presented in the next 
section deals with the finite beam effects neglecting the atmospheric 
effects. The well-known results for the case of infinite beam are derived 
from the finite beam formulation. The latter section develops the t.ieory 
with the effects of the atmospheric turbulence and beam t. uncation. The 
results are expressed as finite integrals. These integrals require numer- 
ical integration. The integrals were programmed on the UNIVAC 1108. The 
program ran into unexpected computer memory problems and the results are 
fragmentary at this time. 

Another objective was r o gainer data on the refractive index structure 
constant using the CA-9 Space Averaging Anemometer. It took Feveral 
months to obtain clearance from NASA to purchase the Anemometer and the 
associated equipment. Campbell Scientific, Inc. took almost 3 months to 
deliver the system. The system is in the process of being set-up over 
the roofs of Patton Hall and Carver Complex on the campus of Alabama A 8 M 
University. 


A paper entitled "Loser Beats Propagation in Turbulent Atmosphere" is 
at the stage of galley proofs and will be published in the forthcoming 
issue of Proceedings C, Indian Academy of Sciences, Bangalore, India. 
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FINITE BEAM THEORY OF ATMOSPHERIC LDV SYSTEMS 


lletcrcdyne detection of optical radiation scattered front aerosols 
has been the subject of several papers. We are concerned specifically 
with tne detection of backscattered radiation and the extraction of 
velocity information of the target from the Doppler shift of the wave 
frequency. The NASA baser Doppler System is of coaxial configuration with 
the same telescope acting as transmitting and receiving antenna. The 
intensity of laser radiation leaving the transmitting telescope has a 
Gaussian profile whose wings arc cut off by the finite size of the tele- 
scope. It is simpler mathematically to treat the beam as an infinite 
Causslan distributed field and Ignore the finite size of the transmitter 
and receiver. Sonnenschcin and Harrigan [1] analyzed the signal-to- noise 
relationships for a coaxial system that heterodynes the backscattered sig- 
nal from atmospheric aerosols assuming infinite Gaussian fields. The real 
systems are of finite size and the radiation becomes a truncated Caussian 
beam. The purpose of this work is to investigate the effect of the finite 
aperture size of the coaxial system on the signal-to-noise ratio and deter- 
mine the conditions when the infinite wave analysis is valid. 

There are a few research papers which deal with the effects of finite 
aperture size in practical systems. We refer to the papers relevant to 
the present work in the following. Classical heterodyne detection of an 
incoming optical signal by superposition of a beam from local oscillator at 
the finite receiver aperture is discussed by Fink [2] and Cohen [3] for plane 
or airy signal fields and by Takonaka, Tanaka and Fukumitsu [A) for Gaussian 
fields. Mandcl and Wolf £ 5 ] considered the problem of detecting a coherent 
light beam in the form of a plane wave or Caussian wave with an aperture of 
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circular or rectangular geometry. Halldorsson and Langerholc [6 J considered 
the lidar set-up and calculated the form factors for atmospheric backscatter- 
ing of laser radiation for misalignment of the axes and axial displacement 
of the aperture. Fried [7] considered the heterodyne detection of atmos- 
pherically distorted wavefront by a finite detector and showed that there is 
a limit to the useful collector diameter for an optical heterodyne system. 
Finally , Lutomirski and Buser [b] analyzed the mutual coherence function for 
a finite optical beam which contains the effects of source geometry and 
coherence loss in atmospheric propagation. The analysis of the coaxial system 
with finite aperture neglecting the atmospheric propagation effects is pre- 
sented in the following. 

The laser beam is transmitted along the Z axis and scattered back at 
the range point by the natural aerosols. The transmitter, receiver and the 
diffuse scatterers along with the coordinate system are shown in Fig. 1. The 
transmitter and receiver lens are in the xy plane and the scatterers are at a 



Fig. 1. Sketch of the telescope, scatterers and 
coordinate systems 

We consider the coordinate vector r^ as a two-dimensional vector in the plane 
of the transmitting lens. The vector r^ is in the plane of the receiver lens. 
The vector p is a two-dimensional vector in the xy ' scattering plane at the 
point P. 
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The wavefront leaving the transmitting lens in assumed to have a 
Gaussian amplitude distribution with standard deviation of u and focussed 


at range f from the transmitter. The transmitted wavefunctlon is given by 


vc*,') = u 


_ at f 
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where k ■ 2”/\ is the optical wave number, and r, ■ r , * 

— X 

The wavefunctlon at the range point P before scattering takes place may 
be written using the Huygens-Fresnal principle: 


VC?) j vfO GCV,. V A 
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The term G (p, r,) is the field at the point P due to a point source at r^ 


and is given by 
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for a homogeneous medium where s (p, r^) is the geometric distance between 
— ► • 

the points p and r^« In the paruxial approximation, we may use 
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After subst itut ing for C (p, r^) in eq. (2), we get 
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The Integration in eq (5) is over the area of the receiver. 

The laser bean will be back-scattered at point P. Applying the 
Huygens-Fresnel principle and paraxial approximation to the scattered beam, 
we obtain the wavefunction at the receiver given by 

, >4 - iJLl £ + WP-*,) 2 , i (_<l uut) 

•(£,) V(F, 4 “ 
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where a 1 is the radar crossection of a single aerosol. 

The incident radiation is mixed with the local oscillator signal whose 
wavefunction is given by 



The signal current i B is given by 

^ * lr ) <c; 

•¥ . 

where W (r^) is the aperture function given by 

» i 'S s »/4 
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We will now substitute equations (5), (6) and (7) in eq (8) and obtain the 
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following equation for the signal current: 
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and 


VtM = w V) t 
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The signal power due to scattering from a single particle at P is 
proportional to the square of the current given by eq (10). To obtain the 
total power, we have to integrate over the probing volume. The pulsed 
system senses scattered radiation from a slab of air of thickness AL - ct/2 
in the paraxial approximation where c is the speed of light and T is the 
pulse length. 

* Let I g denote the total current due to scattering from all particles at 

• 2 
a mean distance L. The total power i 6 proportional to I • 
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In order to obtain the total signal power, we have to calculate 1 , which is 

8 

written as a double integral. The ensemble average of the square of the 
signal current is given by: x < C /£ n 2 - / T* ■—* \ 1 5 
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The point r ^ la in tho plane of the transmitter and the point r, is in 
the plane of tho receiver. We will consider the integrals in each of the 
square brackets one after another. We Introduce sum and difference coodln- 
ates and make use of the circular oymmetry of the beam. 

Define 


"Tt, = R*+ Vt 


P* = -/.> 


Substituting for r^ and r^» we obtain the following relations. 

cF-s)*= 


(13a) 
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Using the relations from eq (13) and (14) , it is possible to write 
the following: . _ 
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( 15 ) 









J is Bess3l function of zero order, 
o 

Manipulations similar to those in eq (13), (14) and (15) yield the 
following results for the remaining integrals in eq (12). 

Define sum and difference coordinates: 
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The Inner integral on the right hand aide of eq (18) is the integration of the 


function 


i 
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over the area of overlap of two circlea, each of diameter D, with centera 
displaced at opposite ends of the vector 


Fig. 2 shows the area of overlap and the relevant vectors. 0 and O' 
are the j J 
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Fig. 2. Area of overlap 


centers of the circles and the length 00' is equal to From the triangle 

OAB, we obtain the following. 
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The inner Integral may be written aa follows using eq. (19), (20) and 


(21) over the hatched area in Fig. 2. 
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Extending the integrati:n over the entire area of overlap, we obtain 
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Using the results of cq (15) - (22) in eq (12), the square of the 
signal current becomes: 
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The Intervals In eq (23) require numerical integration for finite 
aperture. However, they can be evaluated in closed form when the upper 
limits of integration are stretched to infinity which corresponds to 
an infinite beam and aperture. The infinite wave case is obtained by 

4 *♦ 

setting D ■» ® in eq (IS) and (18) and taking W(r) - 1 for all values of r. 
The result is 


ire result is . . ^ 
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The integrals are evaluated using the following relation: 
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The result is the following: 
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The aerosols are of varying sizes and wc Integrate eq (25) over the 

size distribution n(r) where r is the aerosol radius. The result is the 

replacement of a’ by a which is given by 

oo 

C* = | 
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The radar crossoctlon O is related to the backscatter coefficient B by 


0 ■ AiB* Substituting these relations, eq (25) becomes 

v o- J\ iWa*/ 


From cq (11), the total current from all particles at a mean distance L 
for pulsed system is given by (dropping the angle brackets) 
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■ Energy per pul-’t- is defined by 

T 
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where is the average power. By the definition of the wavefunction, the 
number of photons in the pulse 1 b given by 
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Using the relation “ hvN^,, we obtain the relation 
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The noise power is proportional to the square of the noise current and 
is given by 

J2LH ** U* (29) 
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Wo maio use of eq (26), (27) and 2d) Co obtain: 

, irf T «|/»gr ._** — 

^ r [ , , J^(i- i-j] 
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This is cha same at. that obtained in Rot. i, if wo note that a* - - R*72. w ht»re 
R is tho aperture radius. 

For a finite beam, eq (27) and (20) should be modified as follows: 
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lining eq (31) and (32) after u little ulgebra, the signal- to-noiso ratio for 
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finite beams is obtained as follows: 
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where C - /(;. C represents the truncation of the Gaussian beam. For 
Instance, if C ■ 1, the beam is truncated at the points where the intensity 
has dropped to e * As the value of c increases, greater portions of the 
wings are included in the beam and if e -*• «*», there is no truncation. 
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PROPAGATION THROUGH ATMOSPHERIC TURLU1ENCE 


The atmospheric turbulence introduces random fluctuations of amplitude 
and phase into the wavefunction. The turbulence effects are usually incorp- 
orated by modifying the term G (p, r) of eq (3) as follows: 
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where V (p, r) describes the effects of the inhomogeneous medium on propaga- 
tion of a r.pherical wave between the points p and r and is given by 


%'L}^) * + * Stfc* ) 


(35) 


X(p, r) represents the perturbation in the log-amplitude of the wave and 
S(p. r) represents the perturbation of the phase between the points p and r 

The Huygens-Fresnel formulation is modified using eq (34) and (35). With- 
out further elaboration, we write the equation for the signal current using 
eq (8) : 
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whereV^ ■ V (p, r^) ■ i represents the perturbations between the points 
r^ and p, and ^ ^ * 2 ^ " ^2 + i S 2 re P resent8 the perturbations between 
the points p and r^. 
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The ensemble average of Che square of Che signal current is given by 
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The ensemble average in eq (37) may be evaluated using the following 


relations : 
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x and y are independent Cause ian-random variables with means of .. and y, and 
a and b are arbitrary constant. 

Using eq (38), the following relation is obtained: 
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C (r) is the 
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□ log-amplitude covariance given by C ( ft ) - \ 

x * y\ ^ ■* 
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with r - | - r 2 | . 

D (r) it. the wave structure function and may be approximated by 

= 2 (x)' ’ 
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The log-amplitude covariance for a spherical wave valid for 6hort and 
long path lengths, and weak to strong turbulence is given by the following: 
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In order to proceed further, a few approximations are made. Using the 
identity 
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the structure functions are approximated as: 
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Th. second approximation 1. to repine. - J |, a „d , . ; () 
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The sum end difference coordinates defined In e, (13) end «3e) will he 
utilized to simplify eq (42). The result is 
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